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Continuum emission around AGB stars at 1.2 mm 
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ABSTRACT 

It is generally acknowledged that the mass loss of Asymptotic Giant Branch (AGB) 
stars undergoes variations on different time scales. We constructed models for the dust 
envelopes for a sample of AGB stars to assess whether mass-loss variations influence 
the spectral energy distribution (SED). To constrain the variability, extra observations 
at millimetre wavelengths (1.2 mm) were acquired. From the analysis of the dust mod- 
els, two indications for the presence of mass-loss variations can be found, being (1) 
a dust temperature at the inner boundary of the dust envelope that is far below the 
dust condensation temperature and (2) an altered density distribution with respect 
to p(r) oc r~ 2 resulting from a constant mass-loss rate. For 5 out of the 18 studied 
sources a two-component model of the envelope is required, consisting of an inner 
region with a constant mass-loss rate and an outer region with a less steep density 
distribution. For one source an outer region with a steeper density distribution was 
found. Moreover, in a search for time variability in our data set at 1.2 mm, we found 
that WX Psc shows a large relative time variation of 34% which might partially be 
caused by variable molecular line emission. 
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1 INTRODUCTION 

When low- and intermediate mass stars enter the AGB phase 
at the end of their lives, the mass-loss rate exceeds the nu- 
clear burning rate and so dominates the subsequent evo- 
lution. In this process, a circumstellar envelope (CSE) of 
gas and dust is formed which expands at a rate of about 
lOkms" 1 (see for example Habing (1996)). 

The AGB mass loss undergoes variations on different 
time scales. The mass loss gradually increases over a time 
scale of several hundred thousands of years to reach a max- 
imum on the tip of the Thermally Pulsing-AGB (Habing 
1996). The pulsations of the central star cause variations 
over several hundreds of days. Strong variations in the 
mass-loss rate, probably related to thermal pulses occur- 
ring every 5-10 10 4 yr (Vassiliadis & Wood 1993), can lead 
to the formation of circumstellar detached shells. Schoier, 
Lindqvist & Olofsson (2005) have recently studied 7 car- 



bon stars with known detached molecular shells. They find 
that the shells are caused by periods of very intense mass 
loss (~ 10 Mq yr ) and that these shells are running into 
a previous low mass-loss-rate AGB wind. After the intense 
mass-losing period, the mass loss decreases to a few 1O _8 M0 , 
on a time scale of a few thousand years. 

The question then arises if these types of mass-loss 
variations can be deduced from the SED. To answer this 
question, we modelled the SED of a carefully selected 
sample of nearby AGB stars, post-AGB stars and super- 
giants. For modelling purposes it proved desirable to obtain 
millimetre observations of the stars in our sample, as there 
are not many observations available in the literature in 
this wavelength range. Still, a considerable part of the dust 
envelope is so cold that it emits its radiation at millimetre 
wavelengths, making these fluxes necessary to constrain the 
models for the SEDs. 



* Based on observations collected at the European Southern Ob- 
servatory, La Silla, Chile within program ESO 70.D-0133, 71. D- 
0072 and 71.D-0600. 

f Scientific researcher of the Fund for Scientific Research, Flan- 
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The outline of this article is as follows. First, the ob- 
tained fluxes at 1.2 mm are presented. Sect. 2 contains in- 
formation on the sample selection, the observations with 
SIMBA and the data reduction with MOPSI. Sect. 3 is de- 
voted to the aperture photometry. In Sect. 4 the constructed 
models for the SEDs are presented and the quality of these 
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models and the indications for variable mass loss are dis- 
cussed. In Sect. 5, the variability at 1.2 mm is analysed on the 
basis of our observations. The conclusions are summarised 
in Sect. 6. 



2 OBSERVATIONS AND DATA REDUCTION 
2.1 Sample selection 

The AGB stars in the sample include M- and C-stars. We 
have also observed some post-AGB stars and supergiants 
that are known to have circumstellar shells. To improve the 
sensitivity to very faint circumstellar emission, the selection 
was restricted to objects with low mm-background emission 
based on the IRAS 60 and 100 |J.m data (typically the CIRR3 
flag, which is the total 100 fxm sky surface brightness, < 
50 MJy sr -1 ) and the sample was further restricted to stars 
within lkpc (see Table 3). There are a few exceptions, such 
as the unique F-type supergiant IRC +10 420. 



2.2 Observations 

The data were taken with the 37-channel hexagonal bolome- 
ter array SIMBA installed at the 15-m SEST telescope, dur- 
ing 3 observation runs: 2002 September 5-10, 2003 May 9-13 
and 2003 July 13-15. The channels have a half power beam 
width (HPBW) of 24arcsec and two adjacent channels are 
separated by 44arcsec. The filter bandpass is centered on 
1.2 mm or 250 GHz and has a full width half max (FWHM) 
of 90 GHz. 

SIMBA works without a wobbling secondary mirror; the 
maps were made with the fast-scanning technique. The scans 
were performed in azimuthal direction by moving the tele- 
scope at a speed of 80 arcsec s~ with a separation in azimuth 
of 8 arcsec. The size of the scans varies between 500 arcsec 
x 360 arcsec and 600 arcsec x 480 arcsec. 4 to 8 consecutive 
scans were made of each source. This procedure was repeated 
several times for each source during the observation runs. 
Although SIMBA has a hexagonal array design, there are 
still gaps in the spatial coverage. Because the telescope is 
set up in horizontal coordinates, the scan direction on the 
sky changes with the hour angle. By assembling scans taken 
at different hour angles, the gaps in the spatial coverage are 
reduced. Of course the assembling of the different scans also 
increases the signal-to-noise ratio. 



2.3 Data reduction 

The MOPSI software was used to reduce the data. Basically 
the procedure described in the section on MOPSI in the 
SIMBA user-manual (Observers handbook SIMBA 2003) 
was followed, but the reduction was further divided into 3 
different steps. In the first step some fundamental opera- 
tions like despiking, opacity correction and sky-noise reduc- 
tion are performed on each of the individual scans. Also, 
the 4 to 8 consecutive scans of each source are put together 
to form one map. In a second step, the maps made during 
different nights are assembled. If the source was detected 
in the individual maps, the maps were recentered before the 
assembling. Out of this map information is gained about the 







Figure 1. The final map of a Sco, resulting from the observations 
of September 2002. The legend is expressed in mjy, the rms-noise 
on the sky background is 27.7 mjy. 



position of the source, which is used in the last step to define 
a proper baseline and to improve the sky-noise reduction. 

For the absolute calibration, scans were made of Uranus, 
typically at the beginning and/or end of each observing 
night. The data reduction is basically the same as for the 
scientific sources, although a few steps were altered because 
Uranus is a relatively strong and extended source. Only low 
order baseline fits were used (with an appropriate base range 
definition) and neither despiking nor sky-noise reduction 
were carried out, as recommended in the manual. The flux 
of Uranus at 250 GHz was estimated with the subprogram 
'Planet' in MOPSI. Table 1 gives the resulting calibration fac- 
tors for the 3 observation runs. As an example, Fig. 1 shows 
the final map of a Sco resulting from the observations of 
September 2002. 



3 APERTURE PHOTOMETRY 

After the data reduction, the fluxes were determined using 
aperture photometry. For each source we searched for the 
'ideal' aperture: the aperture with the highest corresponding 
signal-to-noise ratio. For the small apertures a correction is 
necessary, because one neglects counts from the source that 
were detected outside the aperture. To apply this correction, 
the point spread function or PSF had to be determined. 

The PSF was calculated using the 3 final maps of 
Uranus from the 3 observation runs. We used Uranus and 
not one of the scientific targets because it is such a strong 
source and its angular diameter is well known. Because it 
is an extended source, the maps do not just represent the 
PSF, but a convolution of the PSF with a uniform disk with 
diameter equal to the angular diameter of Uranus. As we 
wanted to model the PSF with a gaussian profile, we fit- 
ted a convolution of a uniform disk with an elliptic gaussian 
profile to each of the 3 maps and determined the best fit 
with a least-squares method. The variances in the x- and 
y-direction of the 3 resulting gaussian profiles are given in 
Table 1, together with the corresponding FWHM (which is 
2\/2 ln2 a) in the 2 directions. The given uncertainties are la 
errors. As can be seen from Table 1, there is little difference 
between the variances of the 3 gaussian fits. The variances 
of the average PSF were calculated as weighted averages of 
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Table 1. The first row gives the calibration factors for each of the observation runs. The rest of the table contains the 
variances of the gaussian profiles for the 3 observation runs and of the averaged PSF, together with the corresponding 
FWHM. 





Scp02 


May03 


Jul03 


PSF 


calibration factor (to mjy) 


62.0 ±2.6 


64.5 ±2.7 


58.6 ±2.4 






10.44 ±0.47 


10.30 ±0.45 


10.49 ±0.73 


10.40 ±0.30 


<Jy 


10.55 ±0.48 


9.92 ± 0.44 


11.13 ±0.78 


10.43 ± 0.31 


FWHMj 


24.7 ± 1.2 arcsec 


24.3 ±1.1 arcsec 


24.7 ± 1.8 arcsec 


24.49 ± 0.71 arcsec 


FWHM,, 


24.8 ± 1.2 arcsec 


23.4 ±1.1 arcsec 


26.2 ±1.9 arcsec 


24.56 ±0.73 arcsec 



Table 2. The scale factors to correct for count loss within small 
apertures. 



aperture 


scale factor 


40 arcsec 


0.99946 ± 0.00014 


35 arcsec 


0.99658 ± 0.00068 


30 arcsec 


0.9841 ± 0.0023 


25 arcsec 


0.9438 ± 0.0054 


20 arcsec 


0.8419 ± 0.0095 


15 arcsec 


0.647 ±0.012 



a x and G y from the 3 fits, which are also given in Table 1. 
The corresponding FWHM in the x- and y-direction of the 
average PSF are in good agreement with the expected value 
of 24 arcsec, which is the HPBW of one channel. 

With the help of the average PSF, the scale factors 
to correct for count losses with small apertures could be 
determined (see Table 2). An aperture of 45 arcsec already 
contains 99.99% of the flux, so we put the scale factors for 
apertures ^ 45 arcsec equal to 1. The scale factors from the 
3 individual gaussian fits were also calculated to determine 
the rms-errors on the scale factors from the average PSF. 

Once these scale factors were known, we continued by 
determining the ideal aperture and the corresponding flux 
for each of the scientific targets. In the user manual of MOPSI 
one can find the commands to determine the number of 
counts within a given aperture, as well as the error. We 
have discovered that the numbers that MOPSI returns are 
not what we understand to be the number of counts and the 
error. MOPSI considers the number of positive counts in the 
integration area as the total number of counts and the num- 
ber of negative counts as the error. Hence MOPSI adds the 
positive noise to the total number of counts, but neglects 
the negative noise. In our opinion it is better to take the 
positive counts minus the negative counts as the total num- 
ber of counts, because after the data reduction one expects 
the noise to fluctuate around zero. 

To calculate the error on this total number of counts, 
we assumed that the noise follows a Poisson distribution, so 
that the error equals the square root of the total number 
of counts. When the ideal aperture is less than 45 arcsec, a 
scale factor had to be applied to derive the total number: in 
that case the error on the scale factor propagates into the 
error on the counts. The total number of counts within the 
ideal aperture also has to be divided by the area of the PSF, 
ft, being 

FWHM,, FWHM,, n 
4 ln2 



The values of FWHM, and FWHM S from the average 
PSF from Table 1 were used. The error on Q, resulting from 
the errors on FWHM, and FWHM,,, was also accounted for 
in the error on the total number of counts. Finally, we also 
took the error on the calibration factors into account (see 
Table 1). 

To make sure that we do not exclude any extended emis- 
sion by using the ideal aperture, we have fitted gaussian 
profiles to each of the detected sources. We have taken ideal 
apertures that contain at least 99 % of the volume under- 
neath those gaussian fits. Besides the gaussian flux profiles 
of the central sources, no other structures are visible in the 
maps, so there is no danger of contamination of the flux by 
other objects, not even for the largest apertures. 

The resulting flux and the ideal aperture for each of the 
science targets are given in Table 3. 



4 MODELLING OF THE SED 
4.1 Methodology 

For all 23 sources that were detected with SIMBA, mod- 
els were constructed for the SED with dusty (Ivezic et al. 
1999). dusty is a code that solves the problem of radia- 
tion transport in a dusty environment. Its input consists of 
four categories of parameters: information about the exter- 
nal radiation source, about the dust properties, the density 
distribution and the overall optical depth of the dust en- 
velope. From this information, the code calculates the dust 
temperature distribution and the radiation field, by solving 
a self-consistent equation for the radiative energy density, 
including dust scattering, absorption and emission. 

For the modelling of our sources, a spherical model was 
used, in which the external radiation comes from a point 
source at the centre of the density distribution. The spec- 
tral shape of the external radiation is determined by a black 
body of a given temperature. For the oxygen-rich dust en- 
velopes, we used silicates to model the chemical composition 
of the CSE. dusty contains the optical properties of 3 types 
of silicates: hot and cold silicates which are discussed by Os- 
senkopf, Henning & Mathis (1992) and a silicate with optical 
constants from Draine & Lee (1984) (further abbreviated by 
SilOw, SilOc and SilDL). For the oxygen-rich envelopes, only 
models with one type of grain were tested, never a mixture. 
For the carbon-rich envelopes in our sample, a mixture with 
varying ratios of amorphous carbon (AMC, with optical con- 
stants from Hanner (1988)) and SiC (with optical constants 
from Pegourie (1988)) was used. The size of the grains was 
fixed at 0.1 [im. The inner radius of the dust envelope is fixed 
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Table 3. General information about the science targets, together with the derived flux and the corresponding ideal aperture. In case 
no detection was made, three times the rms on the sky background is given as an upper limit. A question mark in the column of the 
distances means that no reliable distance estimate was available. References are given below. 



Name 




RA 


Dec 




spectral type 


distance 
(kpc) 


observation 
run 


flux at 1.2 mm 
(mjy) 


ideal aperture 
(arcsec) 


rms 
(mjy) 


AFfiT, 1 922 


17 


: 04 : 


55.1 


-24 


: 40 : 


: 41 


C a 


0.88 e 


Scp02 


93.7 ± 6.2 


50 


8.0 




16 


: 26 : 


20.3 


-26 


: 19 


: 22 


M1 ^Tah-h 6 

ivx _i_ . •jjau u 


O.19/ 


Scp02 


244 4 + 15 1 


35 


9.2 


aSco 


















Jul03 


0Q7 9 + 147 


40 


8.7 


un aaa 7Q 


06 


17 : 


37.0 


-10 


: 36 : 


52 


B8V b 


0.33 9 


Scp02 


QKfi 1 _i_ 91 a 
OJU.l ZLZ Zrl.U 


40 


10.5 


IRAS 15194-5115 


15 : 


: 19 : 


27.0 


-51 


: 15 : 


: 18 


C a 


0.6 h 


Ser>02 


310 6 + 19 

•j \j * \j _i_ J-CiU 


65 


9.8 


o Cet 


02 


: 16 : 


49.1 


-03 


: 12 : 


: 13 


M7IIIc b 


0.11 e 


Ser>02 


I807 + 11 q 

luU. i _1_ ll.o 


55 


8.0 


ill VjcJ_j ouuo 


23 


: 16 : 


42.4 


16 


: 55 : 


■ 10 


C c 


0.95 e 




49Q A _L Q 


90 


9.8 


IK Tau 


03 


: 50 : 


43.7 


11 


: 15 : 


■ 31 


IviUlllC 


0.26 e 


SpnO.9 


IOU.U ZLZ ZrU.O 


115 


8.6 




1 Q 


24 : 


26.7 


11 


: 15 : 


no 


F8Ia & 




otpuz 


Aac\ 7 + 978 


1 nn 


1 1 7 
±1.1 


VY CM -a 

V I 1VL el 


07 


: 20 : 


54.8 


-25 


: 40 : 


■ 12 


M3 /M4TT& 

1VJ.O / lV.L1.LJi 


0.66 e 


Scp02 


uoy. j zlz ij_.Zr 


60 


10.5 


WX Psc 


01 


: 03 : 


48.1 


12 


: 19 : 


■ 51 


M9 b 


0.54 e 


Scp02 


Q90 q _i_ on n 

O^O.O ZLZ ZrU.U 


115 


6.8 


WX Psc 


















Jul03 


1 co 1 _l inn 

±00. J. ZLZ J-U.U 


75 


8.7 


vi ^nn Ani 

V lOUU .ri_(_]l 


20 


: 07 : 


47.4 


-06 


: 25 : 


■ 11 


M b 


0.66 e 


ivi c\i y yjrj 


1 HQ A 4- £{ Q 

iuo.1 zlz u.y 


70 


7.0 


lAiiO UOUUU -ZrOOU 


08 


: 00 : 


32.6 


-23 


: 56 : 


■ 16 


F5e fe 


? 




< 27.7 




9.2 


HF) ^fil 9fi 


07 


: 13 : 


25.3 


10 


: 05 : 


09 


F5Iab b 


2.3^ 


Sor»09 


< 48.0 




16.0 


HD ^61 26 


















Mav03 
iviciy KJ'j 


< 51.8 




17.3 


S Set 


18 


47 : 


37.2 


-07 


: 58 : 


00 


P6 4CN3l d 


0.62 e 


Scp02 


< 120.0 




40.0 


S Set 


















May 03 


< 32.6 




10.9 


HD 1 87885 

J--L-L-' J-O 1 UUJ 


19 


: 50 : 


00.8 


-17 


: 09 : 


: 38 


F2/F3Tah 6 


? 


Scp02 


< 29.5 




9.8 


HD 1 87885 




















< 23.0 




7.7 


HD 1 79821 


19 : 


: 11 : 


25.0 


00 


: 02 : 


: 19 


G51a, b 


l fc 


Scp02 


< 79.4 




26.5 


Tin 1 7Q^91 


















ivi ay uo 


< 78.7 




26.2 


89 Her 


17 


: 53 : 


24.1 


26 


: 03 : 


23 


F2Ibe b 


0.98/ 


ivi ay UO 


uy.j zlz '-t.o 


50 


9.6 


AFGL 4106 


10 : 


: 21 : 


32.4 


-59 


: 16 : 


: 52 


K0 b 


3.3* 


May 03 


96.0 ± 6.4 


60 


9.6 


IRC +10 216 


09 : 


45 : 


14.9 


13 


: 30 : 


40 


C9,5 d 


0.12 e 


May03 


3764.3 ± 221.4 


70 


12.2 


IRC +20 326 


17 : 


: 29 : 


42.5 


17 


: 47 : 


: 27 


C c 


0.79 e 


May03 


8.9 ± 1.1 


50 


10.9 


IRC +20 370 


18 : 


: 39 : 


42.0 


17 


: 38 : 


: 10 


C 


0.60 e 


May03 


58.9 ± 4.2 


50 


16.0 


RHya 


13 : 


26 : 


58.3 


-23 


: 01 : 


25 


M71lle b 


0.13 e 


May03 


94.7 ± 6.3 


30 


12.8 


VHya 


10 : 


: 49 : 


11.3 


-20 


: 59 : 


: 04 


C6,3c-C7,5e(N6e) d 


0.33 e 


May03 


99.5 ± 6.6 


35 


7.7 


WHya 


16 : 


46 : 


12.2 


-28 


: 07 : 


: 07 


ivmiie 6 


0.12 e 


May03 


280.0 ± 17.2 


45 


7.7 


CPD-56 8032 


17 : 


: 04 : 


47.8 


-56 


: 50 : 


: 56 


WC 6 


1.35 m 


May03 


255.4 ± 15.8 


100 


8.3 


IRAS 14331-6435 


14 : 


: 33 : 


07.9 


-64 


: 35 : 


: 03 


BSIabe 6 


? 


May03 


94.2 ± 6.3 


85 


11.5 


UHya 


10 : 


35 : 


05.0 


-13 


: 07 : 


26 


C6.5,3(N2) d 


0.35 e 


May03 


< 30.7 




10.2 


HR4049 


10 : 


: 15 : 


49.9 


-28 


: 44 : 


: 29 


B9.5Ib-II b 


0.67-f 


May03 


< 28.8 




10.2 


AC Her 


18 : 


: 28 : 


09.0 


21 


: 49 : 


: 53 


F4Ibpvar i > 


0.75" 


May03 


< 107.5 




35.8 


a Her 


17 : 


: 12 : 


21.9 


14 


: 26 : 


: 45 


M5lab b 


0.12 e 


Jul03 


147.5 ± 9.4 


45 


9.3 


RScl 


01 : 


: 24 : 


40.1 


-32 


: 48 : 


: 07 


C6,5ea(Np) d 


0.47 e 


Jul03 


61.2 ± 4.3 


35 


4.1 


TXPsc 


23 : 


43 : 


50.1 


03 


: 12 : 


34 


C7,2(N0) d 


0.23 e 


Jul03 


< 26.2 




8.7 



a = Groenewegen, de Jong & Baas (1993), b = SIMBAD, c = Kwok, Volk & Bidclman (1997), d = Samus ct al. (2004), e = Loup ct al. 
(1993), f = ESA (1997), g = Cohen ct al. (1975), h = Ryde, Schoier & Olofsson (1999), i = Jones et al. (1993), j = Jura, Chen & 
Werner (2000a), k = Josselin & Lebre (2001), 1 = Molstcr et al. (1999), m = de Marco, Barlow & Storey (1997), n = Jura, Chen & 
Werner (2000b) 



by specifying the dust temperature at the inner boundary, 
dusty assumes instantaneous dust formation at this radius. 
When dealing with objects where the expansion of the en- 
velope is driven by radiation pressure on dust grains, dusty 
can compute the density structure in the envelope by solv- 
ing the hydrodynamics equations, coupled to the radiative 
transfer equations. In this case, the only input parameter 
required is the relative extension of the dust shell in com- 
parison to the inner radius: this parameter was fixed at 5000 
(resulting in a dust temperature at the outer boundary of 
the dust shell of about 30 K, to meet the temperature of the 
interstellar medium). In these calculations, dusty assumes 
a luminosity of 10 4 L Q , a gas-to-dust mass ratio of 200 and 
a dust grain bulk density of 3 gem -3 . 



To construct the SEDs, photometric data were col- 
lected from the astronomical database VizieR and from a 
report by van der Hucht (2004, private communication) 
that lists photometric data from the literature between 
12 fim and 6 cm for about half of the sources in our 
sample. To further constrain the model, we made use of 
the IRAS low-resolution spectrometer spectra with wave- 
lengths between 8 and 23 [im (obtained from the website 
http:/ /www.iras.ucalgary.ca/~volk/getlrs_plot.html). To 
estimate the interstellar reddening, a program based on the 
model of Arenou, Grenon & Gomez (1992) was used. Every 
model that was tested, was reddened according to the 
average interstellar extinction law for the Johnson system 
from Savage & Mathis (1979). 
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As a starting point, a model with an effective temper- 
ature of 2500 K and an optical depth at 0.5 /im of 10 was 
taken, with a dust temperature at the inner boundary of 
the dust envelope of 1000 K. For a carbon-rich envelope we 
started out with a chemical composition with 90 % amor- 
phous carbon and 10% SiC, for an oxygen-rich envelope 
SilDL was tried first. Ivezic & Elitzur (1997) have made an 
extensive parameter study from which the influence of each 
input parameter on the shape and scale of the SED becomes 
clear. We have done a parameter study ourselves in a smaller 
parameter space, especially to study the influence of the dif- 
ferent parameters on the predicted LRS spectrum, which 
is not clear from the SEDs presented in Ivezic & Elitzur 
(1997). 

Starting from the model described above, the param- 
eters T e ff, to. 5, Ti and the chemical composition were al- 
tered to improve the fit as follows. First the optical depth 
is adjusted since this parameter has the largest influence 
on the model out of the four mentioned above. It does not 
change the overall shape of the SED, a larger optical depth 
merely shifts the complete SED redwards. Then the chem- 
ical composition is determined. For a carbon-rich envelope, 
the strength of the 11.3 /im feature increases as the percent- 
age of SiC is increased. Concerning the oxygen rich case, the 
type of silicate is determined by the relative strength of the 
9.7 /jm and 18 /j,m features. Since the strength of the features 
is also strongly dependent on the value of To. 5, the silicate 
features can put additional constraint on this parameter as 
well. The effective temperature of AGB sources typically lies 
between 2000 and 3500 K. Since this is a small temperature 
range, the effects on the SED are also relatively small. T e g 
affects the shape of the SED blueward of the peak of the 
energy distribution and only has a minor scale effect red- 
wards of the peak. T e g only plays a role in the SED when 
to. 5 < 1, at larger optical depths the external radiation is 
fully absorbed and has no more influence on the SED. Fi- 
nally, the dust temperature at the inner boundary of the 
envelope, Ti, is adjusted. Ti has no influence blueward of 
the SED peak and the shape of the red part of the SED is 
only affected when Ti becomes smaller than about 500 K. 
The LRS spectrum is also very important in determining 
the correct value of Ti, as the slope of the LRS spectrum 
clearly decreases with decreasing values of Ti . 

The optical depth is varied in steps of the order of 
magnitude of to. 5, e.g. around to. 5 = 50, also models with 
to. 5 = 40 and 60 are tried. T e ff is varied in steps of 500 K, 
Ti in steps of 100 K and the amount of amorphous carbon 
in steps of 1 %. Each model is judged by eye, where equal 
importance is given to the agreement with the SED and the 
LRS spectrum. The values from the default model were only 
altered when this brought along a clear improvement of the 
fit to the observations. 

Two of the sources in our modelling sample, HD 44179 
and CPD-56 8032, are said to be surrounded by dusty disks 
(see respectively Cohen et al. (2004) and de Marco, Barlow 
& Cohen (2002)), while there is strong evidence that the 
dust envelopes from VY CMa and 89 Her are not spherically 
symmetric (see respectively (Josselin et al. 2000) and Wa- 
ters et al. (1993)). For these sources no acceptable model 
could be constructed, which is to be expected as dusty 
assumes spherically symmetric mass loss. For one source, 
IRAS 14331-6435, not enough data were available to con- 



strain a model. For all other 18 stars, the SEDs are pre- 
sented in Appendix A, together with the input parameters 
and the derived parameters from the dusty models. To es- 
timate the errors on the input parameters, we varied one 
parameter while retaining the values from the best model 
for the other parameters, until the fit became unacceptable. 
These errors are clearly only estimates, and not formal errors 
due to interdependencies of the different parameters. 

4.2 Discussion of the SEDs 

4-2.1 Quality of the models 

In general, the models produced by dusty explain the ob- 
servations quite well, certainly when we keep in mind that a 
certain scatter of data around the model is to be expected, 
as our sample stars are variable and the data were taken at 
different epochs. 

There are a few discrepancies that occur in several 
SEDs. For IRAS 15194-5115, RHya, a Her and IRC+10420, 
the observations at visual and/or near infrared wavelengths 
can not be explained by the models. For IRC +10 216, this 
problem was solved by increasing the size of the dust grains, 
but this does not work for all problematic cases. When the 
mass loss is small, the spectrum of the central star still has 
a large influence on the spectrum we observe, especially at 
short wavelengths. For these cases the discrepancy can be 
explained by the fact that the central star was modelled with 
a black body instead of using a stellar spectrum suitable for 
the central star. That the influence of the central star is not 
negligible can also be seen from the LRS spectrum, where 
sometimes the atmospheric absorption by HCN and C2H2 
can be noticed at 14 /im. 

There are also significant deviations between the ob- 
served and modelled submillimetre fluxes. The flux excesses 
that are seen here, will be discussed in Section 4.2.2. The 
excess at centimetre wavelengths is discussed separately in 
Sect. 4.2.3. 

4-2.2 Excess at (sub) millimetre wavelengths ( ~ 100 fim - 
1 cm) 

For almost all of the SEDs, the model provides a good fit 
for fluxes at wavelengths shorter than about 100 /jm, but 
for a few of them the longer wavelength fluxes are not in 
agreement with any of the models that were tested. This 
is the case for AFGL 1922, IRAS 15194-5115, WXPsc, RScl 
and V1300 Aql. In this section, various possible explanations 
for the observed (sub) millimetre excess will be discussed. 

The excess emission at (sub)millimetre wavelengths 
can be due to molecular line emission, which is not included 
in the models. IRC +10 216 is known for its rich molecular 
spectrum. For example, Groenewegen et al. (1997) find an 
on source contribution from molecular lines which amounts 
to 43 % within the passband of the instrument used, which 
lies between 220 and 280 GHz and which has an effective 
frequency of 243 GHz. Nevertheless, the DUSTY-model - 
which does not include molecular contributions - that was 
fitted to wavelengths shorter than 500 fim explains the 
(sub) millimetre data for IRC +10 216 well (see Fig. A4). 
The contribution from molecular lines is too small to 
necessitate a different model. Our SEDs are plotted using 
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a log-log scale in units of Wm~ 2 so that the SEDs show 
the physically relevant energy distribution. We conclude 
that the molecular contribution does not cause a flux 
excess on the scale used and therefore does not influence 
the (sub) millimetre modelling. The molecular contribution 
to the flux will vary from source to source, depending 
on a variety of factors, such as the dust-to-gas ratio, the 
temperature and velocity structure in the envelope, etc. 
However, it is unlikely that these factors will change the 
order of magnitude of the molecular contribution. The 
molecular emission is also proportional to the mass-loss 
rate, which is the best determined parameter of influence. 
Since IRC +10 216 undergoes substantial mass loss, we do 
not expect the molecular emission to cause a noticeable 
excess for any other star in the sample, on the scale used in 
the SED. Therefore, we deem it unlikely that the molecular 
contribution will cause a substantial effect in the parameters 
derived from the SED fitting. 

The alternative explanation for the observed ex- 
cesses is an increased amount of dust radiation in the 
(sub) millimetre range, either (1) because the dust grains 
radiate more readily at these wavelengths than is assumed 
in the model or (2) because of the presence of more cold 
dust than in the model. These effects can be caused by (1) 
the presence of larger dust grains or (2) by an increased 
amount of dust farther away from the star. We have 
explored the first possibility by varying the grain size in 
our models for three test cases: IRAS 15194-5115, RScl and 
WXPsc. Starting from the parameter values of the best 
fitting model with grain size 0.1 /jm, grain sizes equalling 
0.2, 0.5, 1, 5, 10 and 50 /jm were tested. Because larger 
grains not only influence the shape of the SED, but also 
shift the SED as a whole redwards, the optical depth was 
decreased to counteract this effect. For IRAS 15194-5115, 
no model could satisfactory reproduce the SED. For the 
other 2 sources, it was possible to fit the SED, including the 
(sub) millimetre fluxes, but these models fail to reproduce 
the LRS spectrum. As an example: a model with r = 5 
and grain size a = 5 /im fits the SED of WX Psc, but the 
strong features in the LRS spectrum are not reproduced. 
On the other hand, a model with r = 10 and a — 0.75 /im 
reproduces the LRS spectrum, but not the SED. Besides 
single sized dust grains we also tried models with grain size 
distributions n(a) oc a~ q for WXPsc. We used a standard 
power index of 3.5 and a flatter size distribution with power 
index 2, grain sizes between 0.01 and 10 fim and r = 5 or 
10, but the fits are even less satisfactory than for the single 
sized grain models. 

Note that these results imply that a single population 
of large grains or a grain size distribution including large 
grains in a spherical configuration can not be invoked to 
explain the (sub) millimetre excess. It does not exclude 
the contribution of large grains in general. The presence 
of an additional shell of large grains can not be modelled 
with dusty. We have used another radiative transfer code, 
MODUST (Bouwman et al. 2000), to explore the effect of 
such a shell for WXPsc. We have taken the inner shell 
from the model derived from SED fitting as presented by 
Decin et al. (2007) and placed an extra shell with grains 
of size 1 or 10 /im at different positions in the envelope 
(Y = 10,100,1000). Our conclusion is that a shell with 



larger grains can only make a significant contribution to the 
SED at (sub) millimetre wavelengths alone, when the shell is 
placed far from the star. The low dust temperature implies 
that the shell must contain a relatively large amount 
of mass, which was modelled in MODUST by increasing 
the mass-loss rate in the shell. Models with mass-loss 
variations will be presented in the next paragraph and 
further discussed in Sect. 4.2.4. Jura, Turner & Balm (1997) 
have argued that large dust particles formed in a long-lived 
disk around HD 44179 can produce the observed millimetre 
continuum, which is in excess of what simple dust models 
predict. Again, such a non-spherically symmetric configu- 
ration in which a region of enhanced density allows large 
grains to form can not be modelled by dusty. Mauron & 
Huggins (2006) studied the morphology of AGB envelopes 
by imaging the circumstellar dust in scattered light at 
optical wavelengths. For V1300 Aql they found no evidence 
of special structures in the envelope, although the envelope 
itself has an elliptical and not a spherical shape. WXPsc 
seems to have a spherically symmetric extended envelope, 
but a strong axial symmetry close to the star. RScl has a 
spherically symmetric detached shell (see Gonzalez Delgado 
et al. (2003)), hence this shell is probably the cause of 
the (sub)millimetre continuum and there is no need to 
invoke the presence of large grains. Feast, Whitelock & 
Marang (2003) argue that IRAS 15194-5115 has a quasi 
spherical dust shell, with small scale irregularities due to 
the ejection of dust clouds. No information was found on 
the geometry of the CSE for AFGL 1922. Hence on the 
basis of observations of the structure of their envelopes, it 
seems that a nonspherical distribution of large grains is an 
unlikely scenario at least for V1300 Aql and RScl. 

To test whether the observed excess is due to dust far 
from the star, we increase the dust density there in the 
model. The density distribution as calculated from hydro- 
dynamics assumes a constant mass-loss rate and can be ap- 
proximated by a distribution p(r) oc r~ 2 . Besides a calcu- 
lation from hydrodynamics equations, dusty also offers the 
possibility to define the density structure with the use of 
power laws. We tested models with p as a function of the 
relative radius y (scaled to the inner radius of the dust en- 
velope) defined as 

p(y) oc ( y ~-l 1 5 y 5 t 

\ y p 2/i < y «S Y 

with p equal to 1.5, lor 0.5 and j/i equal to 10, 100 or 1000, 
i.e. models with a higher mass loss in the past. 
For AFGL 1922, IRAS 15194-5115, RScl, WXPsc and 
V1300 Aql, such a density distribution improved the fit, both 
for the SED and for the LRS spectrum. The fact that the 
density distribution from the hydrodynamics equations is 
not adequate is an indication for non-constant mass loss. 
This will be further discussed in Section 4.2.4. 

We want to mention IRC +20 326 here as well. For this 
source, the hydrodynamical models overestimated the obser- 
vational fluxes. Models using smaller dust grains (a = 0.01 
or 0.001 /im), could not explain the shape of the SED nor 
the LRS spectrum. The best fit is obtained by using a model 
with a steeper density distribution p oc r -3 in the outer en- 
velope, consistent with a lower mass loss in the past. 

The altered density distribution improves the fit to the 
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data, but it can not explain the feature around 1 mm for 
AFGL 1922 (see Fig. Al). This feature might also appear in 
the SED of IRAS 15194-5115 (Fig. A2), but the fewer obser- 
vations around this wavelength make it difficult to tell. A 
search for the origin of this feature is beyond the scope of 
this article, but molecular emission can be excluded due to 
the size of the excess. 

4-2.3 Excess at centimetre wavelengths 

On the basis of the agreement between model and obser- 
vations at centimetre wavelengths, our sample stars can be 
divided into two groups. For group I, which is the largest 
group, it was possible to construct models which provide a 
good overall fit. WHya is a perfect example of this (see Fig. 
A13). In group II, the fit at optical and IR wavelengths is 
very good, but the models can not explain the observations 
at radio wavelengths, a Her is representative for this group 
(see Fig. A15). For a Her the fit is reasonably good up until 
1.2 mm, the first observations at longer wavelengths are in 
the centimetre range and these fluxes are underestimated 
by the model. Although we are only talking about two 
observations, we do think that this phenomenon is real and 
not due to for example a bad choice of model parameters. 
Firstly, this effect can be seen in several SEDs. We see 
similar excesses for IRC +20 370, WX Psc, a Sco and V Hya. 
For all of these stars there are only one or two observations 
at centimetre wavelengths. For oCet and IRC +10 216 
the excess is less obvious, but in these cases there are 4, 
respectively 8 fluxes that lie systematically higher than the 
model predictions. Secondly, different models over a wide 
parameter range were tested, but non could provide a better 
fit than the one shown. We investigated the influence of 
the different parameters in dusty on the resulting models 
and we found that only the density distribution of the dust 
can make a significant difference in flux at millimetre and 
centimetre wavelengths, without affecting the flux at visual 
and infrared wavelengths. Models with density distributions 
defined by power laws as explained in Section 4.2.2 were 
tested; this sometimes improved the fit, but it could never 
fully explain the observed excess. 

In the literature, there are several articles that com- 
ment on a similar flux excess for evolved stars. In Skin- 
ner & Whitmore (1987), models are constructed for the 
outer layers of the atmosphere and the dust envelope of the 
supergiant a Ori. A composite model of photosphere plus 
chromosphere plus dust emission can explain the flux below 
1 mm wavelength, but to explain the excess at millimetre- 
centimetre wavelengths, yet another component (being a 
partially ionised wind in combination with areas of infalling 
ionised material) has to contribute to the flux. In Knapp et 
al. (1995), radio continuum observations from a sample of 21 
evolved stars with high mass-loss rate and an extended CSE 
are studied, in the hope of detecting new planetary nebulae. 
For 4 stars, excess emission is detected that is too weak to 
be attributed to a compact planetary nebula. The authors 
suggested a partially ionised chromosphere as the source of 
this emission. 

Hence one possible explanation of the excesses we ob- 
serve is the presence of a chromosphere, a Sco, which is an 
early M-type supergiant like a Ori, is known to have a chro- 



mosphere, this is also the case for a Her. The presence of 
chromospheres in supergiants is widely accepted, however, 
it is not yet clear if all AGB stars have chromospheres as 
well. Drake, Linsky & Elitzur (1987) performed a radio- 
continuum survey of cool M- and C-giants in a search for 
chromospheric activity. They concluded that at least in some 
red giants, chromospheres are either absent or they contain 
so little ionised gas that their optical depth at centimetre 
wavelengths is less than unity. oCet, WXPsc, VHya and 
IRC +20 370 all fall in this category of cool M- and C-giants. 
Sometimes these stars are said to have a chromosphere in 
the literature, but without convincing arguments. 

In Sahai, Claussen & Masson (1989) several mecha- 
nisms are discussed to explain the centimetre emission from 
IRC +10 216. They concluded that a simple black body in 
combination with dust emission from amorphous carbon 
is not sufficient to explain the observed fluxes. The most 
probable cause of the excess is an increase in dust emission 
at the relevant wavelengths, which can be justified if long 
carbon molecules (such as PAHs) are present in the dust 
envelope. In Groenewegen (1997), where a model for the 
circumstellar dust shell of IRC +10 216 is presented, the 
dust opacity also had to be increased for wavelengths larger 
than 1 mm. There is free-free emission present as well, this 
emission is negligible for wavelengths below 5 mm, but does 
make a significant contribution to the flux at centimetre 
wavelengths. In this case the free-free emission is not due 
to a chromosphere, but to a small partially ionised region 
with temperatures close to that of the photosphere. 

In summary, the most probable explanation for 
IRC +10 216 is an increase in dust emission at centime- 
tre wavelengths as proposed by Sahai, Claussen & Masson 
(1989), which could be caused by the presence of long carbon 
molecules (such as PAHs) in the CSE. a Sco and a Her are 
known to posses a chromosphere, which could explain the 
excess emission. For o Cet, WX Psc, V Hya and IRC +20 370 
there is no clear evidence for the presence of a chromosphere. 
For V Hya and IRC +20 370, which have a carbon-rich CSE, 
the same mechanism as for IRC +10 216 might explain the 
observed excess. 

4-2.4 Mass-loss variations 

As already mentioned in Section 4.2.1, we propose a 
less steep density distribution than the one calculated 
from hydrodynamics to explain the (sub)millimetre fluxes 
for AFGL 1922, IRAS 15194-5115, RScl, WXPsc and 
V1300 Aql. Hofmann et al. (2001) made an elaborate study 
of WX Psc and they also found that a two component model 
with an inner constant mass loss region (p oc r~ 2 ) and an 
outer region with p oc r~ 1,5 gives a better match to their 
observations than a one component model with a uniform 
outflow. The necessity of a density distribution defined by 
power laws can be seen as an indication for variable mass 
loss. All five power laws that were fitted indicate that the 
mass loss has decreased with time. 

For IRC +20 326, a steeper density distribution p oc r -3 
in the outer envelope was necessary to explain the data. 
This indicates that the mass loss has increased with time. 
There is one source, AFGL 4106, where the outer radius of 
the envelope had to be decreased so that the lowest dust 
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temperature amounts to 42 K instead of 30 K or lower for 
the other sources. This might indicate that the past mass- 
loss rate was so low that we can not detect dust at lower 
temperatures. We were unable to model this source with a 
steeper density distribution in the outer envelope instead of 
the density-cutoff. 

There are several mechanisms that might cause mass- 
loss-rate modulations in AGB stars, each one related to 
quite distinct time scales. Modulations induced by stellar 
pulsations would display a period of a few hundred days. 
Thermal pulses occur only once in ten thousand to hundred 
thousand years. A third timescale of 200 to 1000 years is 
derived from the separation between the nested shells seen 
around IRC +10 216 (Mauron & Huggins 1999) and around 
several (proto-)planetary nebulae. Simis, Icke & Dominik 
(2001) suggest that the interplay between grain nucleation 
and wind dynamics leads to quasi-periodic changes in the 
mass-loss rate on this time scale. Soker (2000) proposes a 
solar-like magnetic activity cycle as a mechanism to pro- 
duce the observed geometry of the CSE surrounding these 
objects. 

From the distance where we find the density profile 
changes, i.e. the place in our model where the break-point 
occurs, we can get a rough idea of the duration of the con- 
stant mass-loss phase. Translating the break-points into time 
scales using a typical outflow velocity of 10 kms~ x , we find 
dynamical ages between 700 — 20 000 yr, with an average of 
3500 yr, taken into account only those sources that show ev- 
idence of mass-loss variability. Of course, if the "constant" 
mass-loss sources are taking into account this average dy- 
namical age will increase substantially since they do not 
trace variability out to ~80 000 yr, i.e., the outer radius in 
our modelling. These time scales exclude stellar pulsations 
as a possible cause of the mass-loss-rate changes detected 
in the SED. The lower limit of the kinematic age is still in 
agreement with the third timescale mentioned above. How- 
ever, the concentric shells seen in IRC +10 216 by Mauron 
& Huggins (1999) that have lead to this timescale, do not 
leave their mark in the SED (see Fig. A4). This leaves ther- 
mal pulses as the most plausible explanation for the observed 
mass-loss-rate modulations. 

We want to stress here the importance of 
(sub) millimetre observations to determine the correct 
density distribution. WXPsc is a clear example. Up to 
100 /jm, a density distribution determined from hydrody- 
namics provides a good fit, but the observations at 400 fim 
and 1.2 mm were clearly underestimated by such a model. 
On the basis of these observations, a power law density dis- 
tribution was determined that provided a better overall fit. 
Measurements at centimetre wavelengths can not be used 
in determining the correct density distribution, because 
of the excesses seen at these wavelengths, as discussed in 
Sect. 4.2.3. 

A second type of variation in mass loss is indicated 
by a low dust temperature at the inner boundary of the 
dust envelope. Under normal circumstances one expects 
dust to form once the dust condensation temperature is 
reached. This dust condensation temperature differs for 
different species of dust, but for silicates and amorphous 
carbon 1000 K can be seen as a typical value. When a dust 
temperature of 150 K is necessary to model the SED, as 



is the case for AFGL4106, this means that the dust that 
dominates the SED has temperatures far below the dust 
condensation temperature, so it is located at very large 
distances from the central star. This implies that the past 
mass-loss rate was far greater then the present mass-loss 
rate. For AFGL4106, which is a post- AGB star, the mass 
loss has stopped and we are dealing with a detached shell. 

All in all, it is not possible to derive detailed informa- 
tion on variations in the mass-loss rate from the SED alone. 
First of all, dusty contains a variety of input parameters 
and the models presented here are probably degenerate. To 
determine the correct model for a given star, as many input 
parameters as possible should be determined in an inde- 
pendent way. With interferometric techniques for example, 
the inner radius of the dust envelope can be determined, 
and thus the temperature at the inner edge. Secondly, the 
physics in dusty can be further improved. In the current 
version there is no way to model two dust shells with each a 
constant, but different mass-loss rate. Of course, this would 
even further increase the level of degeneracy of the models. 
We can conclude that the SED can indicate variations in 
mass-loss rate, but it should be seen as a basis for further 
analysis of the object and not as a final stage of the study. 



5 VARIABILITY AT 1.2 mm 

Several sources were observed during 2 different observation 
runs. Only 2 of those sources, being a Sco and WX Psc, were 
strong enough to produce a flux above the detection limit of 
SIMBA. 

From Table 3 it can be seen that the 2 fluxes of a Sco 
agree with each other within the errors. One does not expect 
to see much variation, as the 2 observations are separated 
by 310 d, while the period of a Sco is 350 days (Percy et al. 
1996). 

WXPsc, a Mira variable with a period in the John- 
son K band of 660 d (Samus et al. 2004), does show a clear 
flux variation. If we take our 2 observations to be coincident 
with minimum and maximum light, we find an amplitude 
of 85.1 mjy and a relative variation (the ratio between the 
amplitude and the average flux at 1.2 mm) of 34%. This 
constitutes a lower limit to the true amplitude at 1.2 mm. 
In Table 4 observations at 1.3 mm are listed that were ob- 
tained with the SEST bolometer between March 1992 and 
March 1994 (Omont, private communication). If we shift 
these measurements to 1.2 mm, correcting the flux with a 
power law that was fitted to our model of WX Psc between 
1.1 and 1.4 mm 

_ ( A[min] \ 4 3 

ail.2mm — F f I \ ^ ) 

we find values of 86.1 and 80.4 mjy. This points towards 
an even larger amplitude and relative variation. In the lit- 
erature, there is not much information to be found on flux 
variations at such long wavelengths. We mention here a re- 
sult on the variability of IRC +10 216, because it is also a 
Mira variable. In Sandell (1994) a figure is shown of the light 
curve at 1.1 mm. We searched for a period in these data us- 
ing Fourier analysis and we found 655.3 d as a period from 
the Lomb-Scargle periodogram. A fit with the following sine 
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Table 4. Fluxes at 1.3 mm from WXPsc and IRC +10 216 that 
were obtained between March 1992 and March 1994 with the 
SEST bolometer (Omont, private communication). No errors were 
derived. 



Name 


flux 
(mjy) 


epoch 
(JD) 


WXPsc 


61 


2448891 




57 


2448941 


IRC +10 216 


2218 


2448695 




2363 


2448764 




2073 


2448942 




1540 


2449070 




1327 


2449164 




2417 


2449315 




1965 


2449440 



function 

S(mJy) = 2329.9 + 362.1 sin f 2tc — "J 447301 ^ 

\ 655 J 

has a variance reduction of 79 %. Here, the relative varia- 
tion is 16 %. From these values we have to conclude that 
the relative variation found for WXPsc at 1.2mm is rather 
large. We also analysed the 1.3 mm data from IRC +10 216 
that are listed in Table 4. We fitted a sine function with the 
period found at 1.1mm, this gives a variance reduction of 
77 % and a relative variation of 25 %, which is about 10 % 
larger than the relative variation at 1.1mm. Molecular line 
emission could explain this difference, as it contributes sig- 
nificantly to the flux at 1.2-1.3 mm: as already mentioned 
in Sect. 4.2.2, Groenewegen et al. (1997) find an on source 
contribution from molecular lines of 43%. If the molecular 
line emission is variable, the observed variations are a su- 
perposition of the variations in the dust emission and in the 
molecular line emission. The larger relative amplitudes at 
1.3 mm could be caused by variability in the molecular line 
emission. As the contribution to the flux from molecular 
lines within the 1.2 mm SIMBA passband is significant, this 
could also explain the large relative amplitude for WXPsc 
at 1.2 mm. 



6 CONCLUSIONS 

AGB mass loss undergoes variations on different time scales. 
In this article, we have addressed the question whether these 
variations also leave their signature in the SED. In order to 
find the answer, models for the SEDs were presented for a 
sample which contained mostly AGB stars. As there are few 
observations available in the millimetre wavelength area, the 
selected targets were observed at 1.2 mm with SIMBA. 

The presented models were constructed with dusty. Al- 
though there are some discrepancies between observations 
and model that occur in several SEDs, there is a good over- 
all agreement between models and observations. Five out 
of 18 sources show a flux excess at (sub)millimetre wave- 
lengths. A modified density distribution compared to a con- 
stant mass-loss rate, expressed in terms of power laws and 
indicative of mass-loss variability, improves the agreement 
for all five sources. This is why we propose that these five 
stars have experienced such strong variations in their mass 



loss that this is detectable in their SED in the form of a 
(sub) millimetre excess. The density distributions that were 
finally used to model these sources, show that the mass-loss 
rate was higher in the past for all five stars. The presence 
of large grains in areas with increased density due to asym- 
metries in the dust envelope can not be ruled out as an 
alternative explanation, although observations of R Scl and 
V1300 Aql show no such asymmetries in the envelopes of 
these stars. One other source, IRC +20326, required a den- 
sity distribution which implies a lower mass-loss rate in the 
past. Besides these altered density distributions, a second 
indication for variations in the mass loss can be found in a 
dust temperature at the inner boundary of the dust envelope 
that lies far below the dust condensation temperature, which 
indicates the presence of a detached shell. We conclude that 
the SED can contain clear indications for variable mass loss, 
but that it is not possible to derive detailed information on 
the basis of the SED alone. This is due to the uncertainties 
on the many free parameters in the dusty models and be- 
cause there are not many options in dusty to adequately 
model a variable mass loss. 

We also analysed the time variability in our 1.2 mm data 
for a Sco and WX Psc, which were observed during two dif- 
ferent observation runs. We found that WX Psc shows a rel- 
ative variation of at least 34%, which is rather large when 
compared to the relative variation of 16% for IRC +10 216 
at 1.1 mm. We propose that a contribution to the flux from 
variable molecular line emission might explain the large rel- 
ative variation of WXPsc at 1.2mm. 
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fitting models. The caption of Fig. Al contains the expla- 
nation of the symbols used, for all other figures only the 
model parameters are given. We comment here on some dis- 
crepancies between LRS spectrum and model that have not 
already been discussed in the article. Several LRS spectra 
show signs of the atmospheric absorption by HCN and C2H2 
at 14 urn. This is the case for IRAS 15194-5115, AFGL3068, 
IRC +10 216 and RScl. At 8^m the signature of the SiO 
fundamental band is visible in the spectra of oCet, IKTau, 
RHya, WHya and aSco. In the spectrum of AFGL3068, 
SiC goes into self absorption around 12 /im, which we should 
be able to model with a larger optical depth. However, a 
larger value of To. 5 degrades the quality of the fit to the 
SED and to the wings of the LRS. This apparent contra- 
diction could indicate a non-spherical mass loss. The dis- 
crepancies between LRS spectrum and model for RScl are 
mainly caused by molecular absorption by HCN at 8 fim and 
around 14/u,m in combination with C2H2. The spectrum of 
RHya shows a strong molecular band from H2O to the left of 
the silicate feature at 9.7 /im. And the discrepancy between 
the model and the blue wing of the LRS for AFGL 4106 and 
IRC +10 420 is probably caused by a lack of Fe in the model. 



APPENDIX A: DUST MODELS 

Table Al gives an overview of the input parameters used 
for the models, Table A2 does the same for the parame- 
ters that were derived from the models. Figures Al to A18 
show the SEDs and the LRS spectra together with the best 
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Table Al. This table gives an overview of the input parameters used for the models. A v is the interstellar reddening derived 
from the model of Arenou, Grenon & Gomez (1992), to. 5 is the optical depth at 0.5 /im, T e g is the effective temperature of the 
central black body and T\ is the dust temperature at the inner boundary of the dust envelope. When no error is given for T e g 
or Ti, this means that the temperature is undetermined. The oxygen rich dust composition is undetermined when to. 5 ^ 0.1, 
except in the case of a Sco. The last column lists any other parameters that were altered with respect to the default model (see 
Sect. 4.1): density structure with parameters p and j/i, relative outer radius of the envelope Y or grain size a. For Y, only the 
order of magnitude is determined. 
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Table A2. This is an overview of the output parameters of the presented models. L 
is the luminosity, n is the distance between the dust envelope and the central star, 
ri/Vc is the ratio of this distance to the radius of the central star. The hydrodynamics 
calculations give a value for the mass-loss rate M and the terminal outflow velocity 
V e . If the terminal outflow velocity falls below 5kms~ 1 , the values of M and V e 
can not be trusted (indicated with an exclamation mark), if T e ff drops below a 
critical temperature inherent to the model, all values listed in the table are uncertain 
(indicated with an asterisk). 
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Figure Al. Ti = 1200 K, AMC = 93 %, t . 5 = 50, T eS = 2500 K, 
p(y) oc y~ 2 between Y = 1 and Y = 100, p(j/) oc y -1 ' 5 between 
Y = 100 and Y = 5000. The upper figure shows the photometric 
data (asterisk) and the model (full line). Error-bars are shown, 
but most of the time they fall within the symbols. A reversed 
triangle represents an upper limit. The lower figure shows the LRS 
spectrum (plus signs) and the model (full line). As a comparison, 
the dotted line shows the corresponding hydrodynamical model. 
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Figure A2. Ti = 900 K, AMC = 96%, r . 5 = 15, T eS = 2500 K, 
p{y) oc y~ 2 between Y = 1 and Y = 100, p(y) oc y~ 1,5 between 
Y = 100 and Y = 5000. 
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Figure A3. Ti = 600 K, AMC = 93%, r . B = 80, T eS = 2500 K. 
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Figure A4. Ti = 600 K, AMC = 97%, r . 5 = 15, T cff = 2500 K, 
a = 0.15 ^m. 
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Figure A5. Ti = 900 K, AMC = 90%, r . 5 = 10, T cff = 2500 K. 
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Figure A6. Ti = 700 K, AMC = 93%, r . 5 = 3, T cB = 2500 K. 
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Figure A7. 7\ = 1000 K, AMC = 90%, r . 5 = 0.5, T cS = 
2500 K, p(j/) oc y- 2 between Y = 1 and Y = 100, p(y) oc j/" 1 
between V = 100 and Y = 5000. 
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Figure A8. Ti = 700 K, AMC = 97%, ro.s = 20, T eS = 2500 K, 
p(y) oc y~ 2 between Y = 1 and Y = 100, p(j/) oc between 
Y = 100 and Y = 5000. 
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Figure A9. Ti = 1000 K, SilOc, r . B = 1, T eH = 3000 K. 
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Figure A10. Ti = 1000 K, SilDL, t . 5 = 20, T off = 2500 K, 
a = 0.15 fim. 
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Figure All. Ti = 1000 K, SilOw, t . 5 = 70, T off = 2500 K, 
p(y) tx y~ 2 between Y = 1 and V = 100, oc j/~ 1,5 between 
y = 100 and Y = 5000. 
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Figure A12. Ti = 300 K, SilDL, r . 5 = 0.1, T cff = 2500 K. 
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Figure A13. Ti = 1000 K, SilDL, r . 5 = 0.1, T eff = 2000 K. 
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Figure A14. Ti = 1000 K, SilOc, t . 5 = 50, T off = 2500 K, 
p{y) oc j/ -2 between Y = 1 and V = 100, p(j/) oc j/ -1 " 8 between 
Y = 100 and Y = 5000. 
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Figure A15. Ti = 1000 K, SilDL, r . 5 = 0.01, T eff = 3000 K. 
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Figure A16. Ti = 150 K, SilOw, t . 5 = 5, T cff = 7000 K, Y 
10. 
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Figure A17. Ti = 1000 K, SilDL, t .s = 0.1, T cB = 3000 K. 
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Figure A18. Ti = 400 K, SilOw, r . 5 



= 10, T eS = 7000 K. 



